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Abstract. We report here theoretical investigations on the complexity of nuclear 
structure, which have been carried out in the framework of different many-body ap- 
proaches, typically applied to nuclear matter and quark matter studies. The variational, 
functional and perturbative scheme are illustrated in their latest developments. The 
effect of various nucleon-nucleon interactions are tested, particularly in the context of 
the nuclear response functions, against a large body of experimental data. The proper- 
ties and decay widths of hypernuclei are shortly revisited, while the equation of state 
of isospin asymmetric nuclear matter leads us toward nuclear systems of astrophysi- 
cal interest, like the neutron stars. Finally the transition from hadronic matter to the 
deconfmed phase of quark-gluon plasma endorses the application of many-body and 
field theoretical techniques to a system with subnucleonic degrees of freedom. 



1 Introduction 

The subject of this report covers a wide range of items and problems in the 
theoretical description of nuclear systems, the common link among them being 
the strong interaction intervening between the constituents of composite systems 
and the many-body complexity of the latter. 

Nuclear and nucleonic structure has been experimentally investigated us- 
ing probes of quite different nature, whose interaction with the nuclear system 
involves forces of different nature, from the electroweak one to the most compli- 
cated strong interactions between relativistic heavy ions. 

We will discuss here at length the information which can be obtained through 
unpolarized and polarized electron-nucleus scattering, briefly touching also neu- 
trino scattering. Strongly interacting probes, like mesons (tt, K) and nucleons 
will be considered in connection with the specific spin-isospin channels they can 
excite in a nuclear system. 

Finally we will briefly illustrate the problems connected with subnucleonic 
degrees of freedom, in particular the transition from nuclear matter to quark 
matter: the so-called dcconfincmcnt phase transition, perhaps the only one oc- 
curring in the early stages of the Universe which one can hope to reproduce in 
the laboratory. 

Different theoretical schemes are employed to test the performances and lim- 
itations of various models for describing nuclear dynamics and spectra. Indeed 
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in the nuclear medium many degrees of freedom come into play and require 
different techniques to be enlightened. 

For example, sticking to the case of electromagnetic probes, one has to con- 
sider not only the one-body nucleonic current, but also the two-body ones, which 
involve pairs of nucleons correlated via the exchange of charged mesons. In the 
domain of deep inelastic scattering the high momentum photons exchanged be- 
tween the electron and the nucleus shed light on the quark structure of the 
nucleon, namely on the complications associated with the subnucleonic degrees 
of freedom and colour confinement. 

Just mentioning all the above items and themes can give the impression of 
many different, disconnected fields of research: on the contrary, theoretical and 
experimental nuclear physics can be viewed as a "magic circle" which involves, 
with various and strong interrelations, all different subfields. One can start, e.g., 
from the single nucleon to proceed toward the structure of nuclei and hypernuclei, 
down to the ideal "nuclear matter" system, then to neutron stars and astrophys- 
ical systems, in which realistic models envisage the existence of a quark-gluon 
plasma phase. And then back, through the confinement mechanism, to nucleons 
and nuclear systems. 

2 Weakly interacting probes: electrons and neutrinos 

The main advantage of electron-nucleus scattering stems from the fairly good 
knowledge of the basic electromagnetic interaction vertex; moreover the virtual 
exchanged photon can penetrate well inside the nuclear system before being 
absorbed, thus involving into the process any nuclear constituent in the whole 
volume of the nucleus. This probe allows to investigate, depending upon the 
four-momentum transfer, both nuclear and nucleonic structure. In the second 
case, which concerns deep inelastic scattering experiments, the partonic structure 
and perturbative QCD corrections can be tested, together with various quark 
models for the confined hadronic states. 

In the former case, which refers to relatively smaller momentum transfers, 
nuclear structure studies can be performed by comparing the theoretical evalu- 
ation of the inelastic electron scattering cross sections with the measured data, 
in a wide range of energy and momentum transfers as well as in different ex- 
perimental setups, appropriate for inclusive (e, e'), one-exclusive (e,e'p), etc. 
measurements. 

In all situations we test the theoretical description of complex nuclei in con- 
nection with both the many-body scheme employed and the role played by the 
NN interaction in shaping cross sections and amplitudes for the excitation of 
specific nuclear states. Taking advantage from the fact that electrons interact 
with the whole nuclear volume, one expects that surface effects are of minor 
importance: thus it is convenient to perform calculations in the framework of 
nuclear matter. 

Both perturbative and variational techniques have been employed to test 
different model descriptions of: 
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1. the two (and three)-body nucleon-nucleon interaction 

2. the nuclear electromagnetic currents and form factors (including the single 
nuclcon current as well as the two-body, meson exchange currents) 

3. the weak nuclear and nucleonic currents (e.g. to investigate the strange form 
factors of the nucleon) 

4. relativistic effects both in the kinematics and in the currents. 

All the above items concur in determining the so-called nuclear response 
functions, which can be directly confronted with the ones extracted from the 
experimental data. 

Different types of response functions can be obtained from the cross sections 
for inelastic electron scattering; we shall list below a few items which have been 
extensively explored. 

Inclusive scattering of unpolarized electrons off nuclei 

In Born approximation the differential cross section (with respect to the 
energy, e', and the scattering angles, Q = (9,<p), of the final electron) is 
given by: 

dQde' ^ a M i v LRL(q,cj) + v T R T (q 7 uj)} (2.1) 

where gm is the Mott cross section, q 2 = u 2 ~q 2 the squared four-momentum 
transfer, v L = (q 2 /q 2 ) 2 , vt = [tan 2 (9/2) — q 2 /2q 2 ] and Rl,t the longitudi- 
nal and transverse (separated) electromagnetic nuclear response functions. 
Relevant physical issues are connected with these response functions, in par- 
ticular the so-called Coulomb sum rule (from the longitudinal response inte- 
grated over the energy) and the scaling properties, both in non-relativistic 
and in relativistic regimes. 
Inclusive scattering of polarized electrons (off unpolarized targets) 

From the difference between the e-nucleus cross sections, with polarization 
of the electron parallel and antiparallel to the incident beam it is possible to 
measure the asymmetry: 

d 2 a+ - d 2 o- 



A = 



d 2 (T+ + d 2 G- 



= Au V lRav(.<1^) + v T R^ v (q,u) + v T >Ry A {q,uj) 
~ ° v L R L (q,uj) +v T R T (q,uj) 

which contains, in the numerator, the parity violating (PV) nuclear response 
functions, stemming from the interference between the nuclear electromag- 
netic and weak (vector and axial) neutral currents. In the above 

G fQ 2 _ 1 4 ( _ Q 2 



A = V- » 3.1 x 1(T 4 t r = -|— (2.3) 

where Gf is the Fermi constant, a the electromagnetic coupling constant, 
Q 2 = -q\ > and v T > = ^tan 2 (9/2) - (ql/q 2 ) tan(0/2). It is worth men- 
tioning that the PV response functions can be sensitive to novel aspects of 
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the nuclear dynamics, different from the ones explored by the electromag- 
netic vertex. In particular they give access to the isoscalar NN force and to 
the neutrons' distribution in nuclei. 
Exclusive versus inclusive electron scattering 

One-exclusive processes, like (e,e',N), in which one nucleon is observed in 
coincidence with the outgoing electron, are more sensitive tests of the nuclear 
models. Indeed the probability of a nucleon to be emitted strongly reflects, 
for example, the momentum distribution inside the nucleus; striking differ- 
ences are obviously found by employing models like a realistic shell model, 
the so-called Hybrid model or even the Fermi Gas model (the latter being, 
manifestly, not appropriate to deal with the details of nuclear structure) . 
A variety of nuclear response functions can be separated by an adequate se- 
lection of the kinematics; in particular some interference between the longi- 
tudinal and transverse components of the electromagnetic current contribute 
to the one-exclusive cross sections, at variance with the inclusive ones. Thus 
one can test different components and matrix elements of the nuclear elec- 
tromagnetic current, including delicate questions like the off-shellness of the 
nucleonic current in the nucleus, the contribution of meson exchange cur- 
rents, the limitations of non-relativistic approaches. [1, 2] 
Finally (e, e', N) is fairly sensitive to the strong interaction of the ejected 
nucleon with the residual nucleus, the so-called final state interaction (FSI), 
which can be responsible for a sizeable distortion of the outgoing nucleon 
wave. Recently a semiclassical approach has been proposed[3] to describe this 
specific aspect: it will be shortly mentioned here, since it provides a natural 
bridge between nuclear matter and finite nuclei. The main research activity 
on the electromagnetic interaction in light nuclei is reported elsewhere [4]. 

Coming back to the main subject of this report, we shall now concentrate 
on the different theoretical methods which have been employed for an accurate 
evaluation of the electromagnetic response functions in nuclear matter. They can 
be roughly grouped into three categories: 

1. The variational method, which is based on the Fermi Hypernctted Chain 
(FHNC) scheme 

2. The functional method, which employs the path-integral approach for de- 
veloping a consistent bosonic loop expansion 

3. The traditional perturbation theory, which encompasses the most frequently 
utilized models, from the Fermi gas (non-relativistic and relativistic) , to the 
Hartree-Fock (HF) model combined together with the important correlations 
provided by the Random Phase Approximation (RPA). 

In concluding this section we briefly mention the case of neutral current 
(NC) neutrino- and antineutrino-nucleus scattering: these processes offer the 
possibility to extract information on the strange form factors of the nucleon. For 
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this purpose two quantities have been considered, the v — v asymmetry [5]: 

/ \NC i \NC 

(a) — er - 

\* I vn — >[i~p \® > vp — >/^+n 

and the ratio of the cross sections for inelastic NC scattering of v(v) on nuclei, 
with the emission of a proton or, respectively, of a neutron: 

T^v{v) _ K^JuiyhP Icy r\ 

K p/n ~ , ,NC ■ [ - Z - b > 

Calculations of both (2.4) and (2.5) have been performed, in an energy range 
from 200 MeV to 1 GeV, within two relativistic independent particle models 
(Fermi gas and shell model); the final state interactions of the ejected nucleon 
have been taken into account through relativistic Optical model potentials. [6, 7, 

8]- 

While the values of the cross sections significantly depend on the nuclear 
model (especially in the lower energy range), the NC/CC neutrino-antineutrino 
asymmetry and the ratio p/n show a rather mild dependence on the model and 
allow one to disentangle different values of the strangeness parameters entering 
into the weak NC axial and vector strange form factors of the nucleon. 

We remind that the Q 2 = limit of the axial strange form factor, g s A , is closely 
related to the problem of the proton spin, a long lasting puzzle raised by the 
measurements, in the deep inelastic polarized lepton scattering, of the polarized 
structure function <?i of the proton. Moreover the information on the strange form 
factors of the nucleon extracted fron neutrino scattering is complementary to the 
similar information, which can be obtained in the above quoted PV polarized 
electron scattering experiments. 



3 Electromagnetic response functions 
3.1 Variational method 

The inclusive transverse response function of symmetric nuclear matter at satu- 
ration density, Rt(q,lo), has been studied by Fabrocini[9] within the Correlated 
Basis Function (CBF) perturbation theory. The main goal of this work is to 
ascertain how Rt is affected by the NN correlations, with a special emphasis on 
the MEC contributions. 

CBF calculations are based upon a set of correlated wave functions 



In >= S 



\n >fg 



(3.6) 
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obtained by applying a symmetrized product of two-body correlation operators, 
f(i,j), to the FG states \n >fg- The following effective structure is adopted for 
/: 

f(i,3)= E / (9) (^)0g\ (3.7) 

9=1,6 

Ojjr 1,6) = (1, °* • °h Sij) ® (l,Ti • Tj) (3.8) 

being the tensor operator; f(i,j) depends upon a set of variational pa- 
rameters, which are fixed by minimizing the expectation value of a realistic, 
non-relativistic Hamiltonian on the correlated ground state. The g.s. energy is 
calculated via the FHNC cluster summation technique, using the correlation cor- 
responding to the Argonne Vu + Urbana VII three-nucleon interaction model 
of ref.[10]. 

The Jastrow (scalar) component of (3.8) accounts for the short range NN 
repulsion, while, among the remaining correlators, the most relevant are the 
spin-isospin and the tensor-isospin ones, which mostly stem from the one-pion 
exchange (OPE) long range part of the potential. 

Configurations up to correlated one particle-one hole (lp-lh) intermediate 
states are considered; the spreading due to the decay of particle (hole) states 
into 2p-lh (2h-lp) states is taken into account via a realistic optical potential 
model. 

The transverse response is given by 

R T (q, u) = j E l< >! 2 S ( w - ( 3 - 9 ) 

n 

where the sum goes over the intermediate excited states \n >, with excitation 
energy ui n ; j(q) is the electromagnetic current operator 

j(q)-j (1) (q)+j (2) (q) (3.io) 

sum of the one-body, magnetic current [j^(q] and the two-body exchange cur- 
rents [j^(q)]- For the latter, the Schiavilla-Pandharipande-Riska model[ll] is 
adopted, which satisfies, by construction, the continuity equation with realis- 
tic V14 Argonne and Urbana potentials. Currents due to intermediate Z\-isobar 
excitations are also included. The sharp energy boundaries of the lp-lh FG 
response are smoothed out by an appropriate folding with a width W(w): 

R T (q,u>) = - { du' R%- lh {q,u/)- Z+'L, w > ( 3 - U ) 

7T J [U) — UJ'Y + [W (W')J Z 

where W(w) = lmW (uj)/M* , W Q being the optical potential and M* the mi- 
cleon effective mass. 

Fig. 1 shows the comparison of the nuclear matter response evaluated in 
ref.[9] and the experimental data on 40 Ca; both the impulse approximation (IA) 
and the full calculation (MEC+IA) are shown. The data are taken from from 
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Fig. 1. Transverse response at q = 300 (a), 380 (b), and 570 (c) MeV/c for 40 Ca and 
nuclear matter. See text. 
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ref.[12] (x), ref.[13] (circles) and ref.[14] (black circles). The global contribution 
of the two-body currents turns out to be positive and provides an enhancement of 
the one-body response, ranging from <~ 20% for the lower momenta (300 MeV/c) 
to <~ 10% for the higher ones (about 500 MeV/c). The tensor-isospin component 
of the correlation is crucial to obtain these results. 

Several variational calculations have been performed in the past for the lon- 
gitudinal response function as well, both in nuclear matter and light nuclei. In a 
recent work Amaro et aZ. [15] have tested a simple approximation to deal with the 
short range correlations (SRC) affecting the lp-lh excited states, which build 
up the nuclear response functions. 

The basic idea of the model (previously tested for the ground state properties) 
consists in truncating the CBF expansion in such a way to retain only those terms 
containing a single Jastrow-type correlation line, 



The expansion in powers of h(r) has the property of conserving the proper 
normalization of the correlated many-body wave function. The nuclear charge 
is conserved as well, at variance with other truncation schemes. The correlation 
employed is the scalar component of a complicated state dependent correlation, 
fixed to minimize the nuclear binding energy in a FHNC calculation with the 
Urbana V14 NN potential. 

Fig. 2 shows the diagrams retained in [15]. The comparison between the calcu- 
lation of the longitudinal response function in the full FHNC and in the proposed 
approximation shows that the difference between the corresponding results does 
not exceed a few parts xl0~ 4 , in a fairly large interval of momentum transfers. 
This proves the reliability of the truncated scheme also for dealing with lp-lh 
excited states; it needs to be further tested, however, for the 2p-2h responses. 

3.2 Functional (path integral) method 

The path integral approach offers a consistent theoretical foundation of the nu- 
clear response functions. Indeed it employs the approximation schemes of QFT, 
namely the perturbative approach and the loop expansion, to provide a fully con- 
sistent criterium of selecting classes of diagrams, still preserving the fundamental 
symmetries and invariances of the choosen model Lagrangian. 

Starting from a Lagrangian which describes interacting nucleons, pions and 
eventually p mesons, all coupled to an external electromagnetic field, one derives 
the nuclear polarization propagator and response functions through derivatives 
of the generating functional 



h(r) = / 2 (r)-l. 



(3.12) 




(3.13) 
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Fig. 2. Diagrams considered in the model of ref.[15]. The dotted lines represent the 
correlation function. The oriented lines represent particle and hole wave functions. 
The black circle indicates an integration point, while the black squares indicates the 
integration point, where the charge operator is acting. 

where ip,^ are the nucleoli fields, <P and A^, respectively, the pion and electro- 
magnetic fields and j^f* the external source of the latter. In the above 

2 

— 777 — 

C = - M)V + h(d^f - -f$ 2 - iphs-nl) ■ <P (3.14) 

is the model Lagrangian (without the p meson) and 

J^ = e^i{l + T 3 )j^ + e[$xd^] 3 (3.15) 
Bp, = e 2 g^+$- (3.16) 

the required couplings of nucleons and pions to the e.m. field. The results ob- 
tained, within this formalism, for the longitudinal and transverse responses in 
the inclusive electron scattering arc shown in the work by Cenni et al. [16] 
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It is worth mentioning that an important connection between the variational 
technique and the path integral approach has been recently explored in ref.[17]. 
These authors suggest a new method, based on the variation, within a path 
integral framework, of a trial Hamiltonian: it turns out that a particular choice 
of the latter corresponds exactly to the use of a Jastrow correlated Ansatz for 
the wave function in the FHNC approach. Thus the new formalism generalizes 
the FHNC and CBF techniques, allowing for their extension to relativistic and 
field theoretical problems. 

3.3 Perturbation theory 

Several approaches based on perturbative schemes have been applied to the 
evaluation of nuclear response functions both for inclusive (parity conserving 
and parity violating) and exclusive processes. Most calculations are based on 
non-relativistic Hamiltonians and currents, but attempts to include relativistic 
effects, especially in the electromagnetic currents and vertices, have been done. 

Many-body perturbative techniques have been developed to treat quasielas- 
tic electron scattering in nuclear matter. The Green's function method is par- 
ticularly suited to express inelastic cross sections through the imaginary part 
of the polarization (particle-hole) propagator; nuclear matter results compare 
fairly well with the experimental data for energy transfers in the region of the 
quasielastic peak and above, providing the A33 resonance and mesonic degrees 
of freedom are incorporated in the nuclear matter polarization propagators. 

The role played by the NN interaction in reshaping the nuclear matter re- 
sponse functions with respect to the pure FG one has been widely investigated, 
by employing both phenomenological effective interactions and G-matrix pa- 
ramctcrizations based on the one-boson exchange potential. Starting from the 
naive Fermi gas description, one can easily obtain the Hartree-Fock (HF) or 
Brueckner-Hartree-Fock (BHF) polarization propagators; the next step is then 
to include RPA correlations, both in the "minimal" version of the so-called ring 
diagrams and in the framework of the fully antisymmetrized RPA. 

Numerical calculations of the HF+RPA response functions can require large 
computing time. Indeed, although in nuclear matter the ring approximation for 
the polarization propagator can be analytically evaluated, as soon as one dresses 
the particle and hole propagators with HF self-energies, even the simple HF p-h 
propagator, IIhf{<1i w), requires numerical integrations when the HF self-energy 
is derived from any realistic NN interaction. 

Moreover the simple algebraic equation for the ring polarization propagator 
turns into an infinite series of complex integrals when the exchange matrix ele- 
ments of the p-h interaction are taken into account in the fully antisymmetrized 
RPA, which is diagrammatically shown in Fig. 3. An approximate, but reliable, 
treatment of these contributions is thus compulsory. 

This has been achieved with the method of the continued fraction (CF) 
expansion [18], which will be briefly sketched in the following. The CF-like ex- 




Fig. 3. Diagrammatic representation of the perturbative expansion for the polarization 
propagator in RPA. The dashed line represents the particle-hole interaction V. 



pansion for the polarization propagator can be written as 



n(!°) 

n RPA = — s (3.17) 

1 - A 



where TT^ ^ is the free (or HF) particle-hole propagator and all symbols appear- 
ing in (3.17) arc functions of both the energy and momentum transfers. The 
CF approach at n-th order reproduces exactly the perturbative series express- 
ing the antisymmetrized RPA at the same order while it approximates higher 
orders. Setting, for example, 

oo 

n RPA = ^2 (3.18) 
the first order in CF provides the following approximation for ZZ"'™) : 



7j(») ~ tj(o) 



(3.19) 



where TI^ = Il^iVII^ + n^ ex is the sum of the direct and exchange first 
order terms. With this approximation the summation in (3.18) is trivial, yielding 

77(0) - g(0) ( 3 20) 

UcF1 - 1 - 77( 1 )/iT(°) 1 - 4V/7(°) - 77( 1 )ex/7 7 (o) • ^'^ u ^ 

The next step extends the approximation to the exact second order term; the 
corresponding RPA propagator reads: 

TjRPA _ £^ (o o-i\ 

CF2 i - nw/nw - {i7( 2 ) //7(°) - [nw /tt(°)] 2 } [ ' 
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~ 1 - 4V77(°) - IlW ex /77(°) - {77( 2 ) e * /77(°) - [nW ex /77(°)] 2 } ' 

where the approximation provided by the first order CF expansion has been 
subtracted off. 

Numerical calculations for the longitudinal and transverse responses, using 
a G-matrix based on the Bonn potential, show that the CF expansion rapidly 
converges, giving a few percent correction between the first and the second order 
CF results. On the contrary, the differences between ring and RPA (in the CF1 
or CF2 approximation) responses appear to be sizeable, both in the longitudinal 
and in the transverse channel, with, perhaps, the exception of the spin-isovector 
channel, as it is illustrated in Fig. 3. We also remind here that a very useful 
(and accurate) approximation, consisting in a "bi-parabolic" parameterization 
of the HF self-energy of the nucleon[19], allows to express analytically IIhf, 
thus considerably reducing the computational time. 

Within the perturbative approach it is also possible to deal with relativistic 
effects, which might become important for energy/momentum transfers larger 
than 0.5 GcV or so. Besides the corrections stemming from the relativistic kine- 
matics and vertices, which have been completely taken into account in the frame- 
work of the Fermi gas model, several investigations, based on a fully relativistic 
Lagrangian with nucleonic and pionic degrees of freedom, both coupled to the 
external electromagnetic field, have been performed. These approaches usually 
include both the single-nucleon e.m. current and the two-body MEC and ac- 
count for relativistic effects within an expansion in terms of suitably "small" 
variables, like the energy/momentum transfers in units of twice the nucleon mass 
(k = q/2M, A = w/2M). A new relativistic expansion for the matrix elements of 
one- and two-body e.m. currents has been recently proposed by Amaro et al.[l] 
in terms of the parameters rji = Pi/M (characteristically of the order of 1/4), 
{pi} being the initial-state nucleon momentum inside the nucleus: clearly this 
formalism will be especially suited in the regime of GeV energies, where, on the 
contrary, n and A are no longer small with respect to 1. 

y— scaling and Coulomb sum rule The concepts of y-scaling and Coulomb 
sum rule have been throughly discussed in the past since they offer important 
tests for the accuracy and reliability of the nuclear model description underlying 
the calculation of the electromagnetic response functions. Both of them require 
a careful investigation in order to disentangle, from the full expression of the 
longitudinal and transverse response functions, those dividing factors such that 
the reduced responses have scaling properties or (for what concerns the longitu- 
dinal response) satisfy the Coulomb sum rule. This task is not trivial since, as 
shown by the Relativistic Fermi Gas (RFG) model, the e.m. form factors to be 
divided by are not simply factorizable, due to the structure of the relativistic 
one-nucleon current. 

In the work by Barbara et aZ.[20] the dividing factors Gl, Gt are constructed 
in the RFG for the longitudinal and transverse response functions in such a way 
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Fig. 4. Nuclear matter responses for fci? = 195 MeV/c at q = 500 MeV/c: the dotted 
line corresponds to the free response, the dashed one to the ring approximation and the 
solid line to the RPA-CF1 response. The kinematics is relativistic and the particle-hole 
interaction is derived from a realistic G-matrix. 



that the reduced responses, so obtained, scale; another dividing factor, Hl, yields 
a (different) reduced response which fulfills the Coulomb sum rule (CSR). The 
relevant point here is that Gl, Gt and Hl are all found to be only weakly 
model-dependent, thus providing essentially universal dividing factors. It be- 
comes thus possible to test whether specific models scale and satisfy the CSR. 
This investigation has been carried out for the Hybrid model [2] and the Quantum 
Hadrodynamic modcl[21] (QHM), each one with different types of interaction ef- 
fects which go beyond the strict RFC It is found that while the Hybrid model 
scales and obeys the CSR, the QHM does neither. The Coulomb and higher 
order sum rules have also been investigated by Amore et a?. [22] in the RFG, by 
employing two different methods, namely by exploiting the scaling properties 
of the longitudinal response function and by enforcing the completeness of the 
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states in the space-like domain via the Foldy-Wouthuyscn transformation. 

3.4 Exclusive versus inclusive electron scattering scattering 

The connections between exclusive and inclusive electron-nucleus scattering 
within the framework of the plane-wave impulse approximation (PWIA) have 
been investigated in a work by Cenni et al.[2]. These authors test the interplay 
between the (model independent) kinematical constraint and the (model depen- 
dent) features of the spectral function in providing the exclusive (and inclusive) 
nuclear responses. The RFG and the Hybrid model are employed to provide 
a link between finite and infinite Fermi systems and to assess the impact of 
the confinement of the struck nucleons on the inclusive charge response. One of 
the important outcomes of this work is that, with an energy-shift and rescaled 
Fermi momentum, an effective RFG response can be obtained whose first three 
energy-weighted moments agree quite well with the analogous quantities evalu- 
ated within the (confined) Hybrid model. 

A more recent work[3] evaluates semi-classically the exclusive (e, e',p) cross 
sections, yet partially preserving the simplicity of nuclear matter calculation, but 
going beyond the PWIA. Indeed this approach allows to include the distortion 
of the outgoing nucleon wave (the so-called final state interaction, FSI) at least 
within the mean field approximation. 



Sticking to the IA, the cross section for the (e, e',p) scattering process, which 
is schematically illustrated in Fig. 5, reads: 




Fig. 5. 




(3.22) 
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In the above: 



da \ 1 mi „„ M 2 , N , 

df2 e J eN cos 2 (0/2) ee' E(p)E N 

is the single nucleon cross section, while 

1 M < A\a\6 (q - H - ij) a g \A > 

S («I=(5?1W ( " 4) 

is the nuclear spectral function, which in (3.22) depends upon the missing mo- 
mentum p m = q p N and the missing energy £ = En — u) (see Fig. 5 for the 
other symbols). 

Different nuclear models generally provide different spectral functions, as one 
can see, for example, by comparing the ones obtained in the RFG, 

S RFG (p, £) = j^yO{k F - p)S j£ - yjkjr + M 2 - V^Tm^| (3.25) 

and in the Hybrid model (harmonic oscillator bound states combined with a 
continuum of unbound states): 

S HM (p, S)=^25{S- (N max - N)cj } n N (p), (3.26) 



N=0 



where 



n "(P)= E ^I^WI 2 (3-27) 

n,e=N-2n 



is the momentum distribution for the N-th shell. The spectral function (3.26) 
is illustrated in Fig. 6, where a dashed line is drawn in correspondence with the 
support for the spectral function of the RFG, eq. (3.25), which is infinite along 
that line and zero elsewhere. 

Within the semiclassical framework, the Wigner transform of the spectral 
function reads: 



[S(£)] w (R,p) = e 



2 

e F - -4 V(R) 

2M* v ' 



5{£ 



2 

e F - -4 V(R) 

2M* K ' 



(3.28) 



where V(R) is a suitable shell model potential (e.g. the harmonic oscillator or the 
Woods-Saxon potentials) and e F = k F (R)/2M* + V(R) the local Fermi energy. 
It is shown in rcf. [3] that in the semiclassical approach it is possible to account 
for the distortion of the ejected nucleon (illustrated by the black dot in Fig. 5) 
in a relatively simple way: an appropriate distortion operator is introduced, for 
which two heuristic assumptions arc made, corresponding to the eikonal and 
to the uniform approximations, respectively. The corresponding exclusive cross 
sections turn out to be quite different, hence stressing the relevance of the FSI 
in the analysis of the exclusive processes. 
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Fig. 6. The spectral function in the Hybrid Model. 



4 Strange probes and strange matter 
4.1 JST+ nucleus scattering 

The quasielastic if + -nucleus scattering has been thoroughly investigated by De 
Pace et al.[23, 24, 25]. This subject shares several motivations with the items in 
the previous sections. The first one relies on the relatively small (as compared 
with strong interaction processes) K + N cross sections: as a consequence the K + 
projectiles can penetrate deeply inside the nucleus, similarly to electrons and 
photons, thus allowing to investigate collective effects in the nuclear response. 
The second reason of interest is the possibility of exploring scalar-isoscalar ex- 
cited states, which are the dominant channel for the if + -nucleon coupling. This 
channel can be partly explored in the longitudinal electron scattering response 
functions, but an isospin separation in addition to the customary longitudinal- 
transverse one has never been, till now, carried out in (e, e',X) experiments. 1 
Finally the observed excess in elastic if + -nucleus scattering cross sections seems 

1 Such a separation should become possible with PV electron scattering experiments 
and it will be of extreme interest to compare the latter with the K + quasielastic 
scattering. 
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to require, to be explained, some enhancement induced by the nuclear medium 
itself. 

A simple phenomenological approach, which expresses the above mentioned 
cross sections by means of an effective number of participant (N e ff), 

d2(J AT d<J K+N Dl \ ori N 

dn^ = Neff ^T R{q ^ (429) 

seems to require values of N e ff too large with respect to what expected in 
the Glauber theory; moreover (4.29) fails in describing the observed collective 
phenomena at low momentum transfer. 

The model proposed in [23] allows to evaluate the nuclear response R(q, u>) 
within a continuum RPA, using an effective G-matrix particle-hole interaction 
and a Woods Saxon mean field, complemented by an appropriate spreading width 
of the ph states. The kinematics is relativistic. An improved Glauber theory is 
used for the reaction mechanism, including one- and two-step contributions. 




50 100 150 50 100 150 200 250 
cj (MeV) co (MeV) 

Fig. 7. The K+- 40 C& cross sections at q = 290 MeV/c (left) and q = 480 MeV/c 
(right); data are taken from ref.[26]. For the explanation of the theoretical curves see 
text. 



In Fig. 7 we show an example of the results obtained in this framework: 
the experimental data for K + — 40 Ca scattering are compared with the 1 step 
(free)+ 2 step calculation (dotted line), the 1 step (using the N e ff formula) + 2 
step (dot-dashed line) and the 1 step (RPA)+ 2 step (dashed line). While the free 
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response is clearly inadequate to reproduce the data, the other two calculations 
fail in reproducing the excess cross section appearing in the low energy region at 
small momentum transfers. However if one empirically reduces by about 50% the 
effective ph interaction in the scalar isoscalar channel, then the 1 step (RPA)+ 2 
step turns out to be in very good agreement with the experimental data, in the 
whole range of momentum transfers and nuclei ( 12 C and 40 Ca) explored. This is 
a new, challenging information on the ph force acting in this channel. 

4.2 Hypernuclei (structure and decay) 

In the previous paragraph we have considered the interaction of a strange parti- 
cle, the K + (su), with the nuclear medium. Should a K~(su) hit a nucleus, most 
of the times it gets absorbed by a nucleon, leading to the formation of a hyper- 
nucleus, namely a nucleus containing one (or more) hypcron (sqq). The lightest 
and more easily produced hyperon is the A particle and a great interest is fo- 
cussed on the investigation of the ground state properties and decay mechanisms 
of /1-hypernuclei. Several experimental data are already available on light and 
medium hypernuclei, but a reacher set of measurements are expected from the 
FINUDA experiment, which is just starting operating in the Frascati National 
Laboratory of INFN. 

The weak decay width of /1-hypernuclei can be evaluated with the Green's 
function (polarization propagator) method, which is conveniently implemented 
in nuclear matter and extended to finite systems through the so-called Local 
Density Approximation (LDA) . The width of the hypernucleus is derived from 
the A self-energy 

r A (k,p) = -2hnS A {k,p) (4.30) 

after integrating over the A momentum distribution. 

Three different mechanisms contribute to the weak decay of a A in the nu- 
cleus: 

— the mesonic decay, 

a^itN (r M ), 

which also occurs in free space and is relevant only in light hypernuclei, 
since the Pauli blocking prevents the low-momentum outgoing nucleon to 
be emitted in a medium-heavy nucleus; 

— the two-body non-mesonic (NM) decay 

AN — > NN (A), 

— the three-body non-mesonic decay 



ANN — ► NNN (r 2 ). 
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A 




(a) (b) (c) (d) 




(e) (f) (g) 

Fig. 8. Lowest order terms for the A self-energy in nuclear matter. 

A new evaluation of the A self-energy has been carried out [27], which includes 
the lp-lh and 2p-2h polarization propagators (related to i~i and i~2, respectively) 
in a RPA scheme. Fig. 8 illustrates a few examples of the diagrams whose imag- 
inary part is contributing to the hypernuclear decay width. The particle-hole 
interaction embodies the exchange of tt and p mesons, together with short range 
repulsive correlations described by the two Landau-Migdal parameters, g' and 
g' A (the latter being used when one vertex of the exchanged pion coincides with 
the weak AwN vertex). The nuclear matter calculation is adapted to finite nuclei 
using the LDA, which amounts to replace the constant Fermi momentum with 
a local one, expressed by kp(r) — [37r 2 p(r)/2] 1 / 3 , p(r) being the nuclear density. 
The width associated to the decay of a A with momentum k is then 

r A (k) = J dr\^ A (r)\ 2 r A [k,p(r)} (4.31) 

where ip A (r) is the A wave function in the nucleus. The results one obtains turn 
out to be very sensitive to the latter, the best choice corresponding to a Woods 
Saxon potential well which reproduces the measured s and p energy levels of 
the A in the nucleus. Finally the total width of the hypernucleus is obtained by 
integrating (4.31) over the momentum distribution of the A itself: 



r A = / dk|A4k)| 2 r4k). 



(4.32) 
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The results obtained within this framework are illustrated in Fig. 9, where the 
decay width of several hypernuclei (from the "light" 12 C to the heaviest 209 Bi 
and 238 U) is shown, in units of the free A decay width, /free, as a function of the 
mass number. The separate contributions, r M , I\ and r 2 , are also displayed. 




A+l 

Fig. 9. A decay widths in finite nuclei as a function of the mass number A; the labels of 
the various curves correspond to: the mesonic width (-Tm), the one- (A) and two-body 
(A) induced decay widths, the total non-mesonic width (-Tjvm = A + A) and the 
total (i~r), sum of all independent contributions. 

The theoretical results are in good agreement with the data over the whole 
hypernuclear mass range explored. The saturation property of the AN ~ * NN 
interaction in nuclei clearly appears in the flattening of the decay width as the 
mass number increases. In spite of these satisfactory results, there remains to be 
explained the large experimental value of the ratio r n /r pi which does not seem 
to be in agreement with the present theoretical schemes. 

Concerning the ground state properties of hypernuclei we should men- 
tion here a research project by Co et al[28], who intend to investigate the struc- 
ture of hypernuclei within the FHNC scheme. Starting from a suitable basis of 



Nuclear matter, nuclear and subnuclear degrees of freedom 21 

single particle wave functions (Woods-Saxon and/or HF wave functions with 
finite range interaction) a Jastrow correlated many-body state will be obtained 
by developing a FHNC formalism for a system of A nucleons plus one impurity 
(the hyperon). Binding energy and first excited states can also be calculated. 



5 Asymmetric nuclear matter 

5.1 Equation of state of asymmetric matter 

Nuclear matter is usually considered as an isospin symmetric system, Z = N = 
A/2, where the Coulomb interaction between protons is switched off. In reality 
heavier nuclei are more and more isospin asymmetric, the larger number of neu- 
trons being there to balance the increasing Coulomb repulsion. An interesting, 
ideal extension of these systems is the asymmetric nuclear matter, where the 
neutron excess is measured by the asymmetry parameter / = (N — Z)/A. 

Several investigations have been performed on the ground state properties of 
(strongly) asymmetric nuclear matter, with large asymmetry parameter, up to 
values as large as 0.8. In the extreme limit / — ► 1 one obtains neutron matter, 
which is obviously unbound as far as nuclear forces are concerned, but deserves 
very interesting applications in the description of neutron stars, which are bound 
by gravitational force, but can be considered as a unique realization of this ideal 
system. 

The equation of state of asymmetric nuclear matter contains a symmetry 
term which depends upon the asymmetry parameter and the density of the sys- 
tem; the equilibrium conditions are crucially determined by the NN interaction: 
various effective forces can produce quite different behaviours with p in the po- 
tential energy contribution to the symmetry energy. [29] 

A study of non-equilibrium properties of asymmetric nuclear matter has been 
carried out by Di Toro et a/. [30] starting from two Vlasov equations, for proton 
and neutron liquids, coupled through the mean field. The unstable growth of 
density perturbations is influenced by the initial asymmetry: in neutron rich 
nuclear matter the formation of larger fragments, in which isospin symmetry is 
restored, is favoured. 

At low (subsaturation) densities the system develops collective dynamical 
instabilities of isovector nature, the so-called spinodal decomposition: it cor- 
responds to the growth of density pertubations leading to a liquid-gas phase 
separation in which symmetric nuclei coexist with a neutron gas. The instabil- 
ity region (for example in the [T, p) plane, turns out to be reduced in strongly 
asymmetric systems. 

This phenomenon can be tested in heavy ion reactions with radioactive beams 
and could provide important information on large systems of astrophysical in- 
terest. 



22 Wanda M.Alberico 



5.2 Neutron stars: structure and equation of state 

A large amount of information on neutron stars should be available in the next 
few years from the new generation of X- and 7-ray satellites. In view of this, 
neutron stars are the subject of many theoretical studies, aimed to predict their 
structure on the basis of the properties of dense matter. The equation of state 
(EOS) of neutron stars covers a wide density range, from ~ 10 g/cm 3 in the 
surface to several times nuclear matter saturation density (p = p n <~ 2.8 x 
10 14 g/cm 3 ) in the center of the star. 

The core (interior part) of a neutron star is believed to consist of asymmet- 
ric nuclear matter with a consistent fraction of leptons. At ultra-high densities, 
matter might undergo a transition giving rise to other exotic hadronic compo- 
nents, like hyperons, a K~ condensate or a deconfined phase of quark matter. 
This occurrence, in turn, could critically influence the evolution of neutron stars 
and could eventually lead to the formation of black holes. [31] 

Baldo et al. have derived the static properties of non-rotating neutron stars 
in a conventional framework, using the microscopic EOS for asymmetric nuclear 
matter, derived from the Brucckner-Bethe-Goldstone many-body theory. [32] 
They use the Argonne V14 and the Paris two-body interactions, implemented by 
the Urbana model for the three-body force. The latter are incorporated into the 
Brucckner scheme by reducing them to affective two-body, density dependent 
forces. The parameters are adjusted to reproduce the empirical nuclear matter 
saturation point. 

The calculated EOS allows to compute masses and radii as a function of 
the central density no. Assuming that a neutron star is a spherically symmetric 
distribution of mass in hydrostatic equilibrium, and neglecting the effects of rota- 
tions and magnetic field, the equilibrium configurations are obtained by solving 
the Tolman- Oppenheimer-Volkoff (TOV) equations for the total pressure and 
the enclosed mass: 



l + P(r) \ ( i Anr z P{r) 



dP(r) Gm(r)p(r) \ c 2 p(r) J \ c 2 m(r) 
dr = r 2 / 2Gm{r) 

\ rc 2 

dm{r) 



(5.33) 



4irr 2 p{r), (5.34) 
dr 

where G is the gravitational constant. The maximum mass configuration ob- 
tained in ref. [32] range from 1.8 to 2.13 solar masses, while the corresponding 
radii vary from 8 to 10.6 Km. These values are consistent with the present ob- 
servations. 

A challenging suggestion is found in a work of Bombaci[33], concerning the 
semiempirical mass-radius relation extracted for the X-ray burst source 4U 
1820 -30: the latter cannot be reproduced by neutron star models based on con- 
ventional EOS. It is thus suggested that the onset of a phase transition to a K~ 
condensate could explain, with a suitable choice of the parameters, the observed 
mass-radius relation. 
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5.3 Neutron stars: the effects of superfluidity 

As already mentioned in the above subsections, it is believed that the inner 
crust of a neutron star, at a density of about 3 x 10 11 g/cm 3 (corresponding to 
about 10~ 3 p ), is a superfluid and inhomogeneous system, consisting of a lattice 
of nuclei immersed in a sea of neutrons and an approximately uniform sea of 
electrons. Such a configuration persists up to roughly po/2. At low temperatures, 
this system is superfluid with a positive Fermi energy ep\ the estimate of the 
pairing gap in this inhomogeneous medium has been carried out by Barranco et 
al. [34] . They solve the Hartree-Fock-Bogoliubov (HFB) equations in a Wigner- 
Scitz cell, with a Woods-Saxon potential; the gap is calculated self-consistently 
without altering the single particle levels. Use is made of the Argonne two-body 
interaction. 

The HFB method is suitable to describe the halo-nuclei and/or neutron- 
rich nuclei. Surface collective modes (phonons) mediate the induced interaction 
producing Cooper pairs and a bosonic condensate. 

In a second work[35] the quantum calculation is compared to the LDA. Here 
the two-body interaction is assumed to be a Gogny force and is included in the 
HFB equations via the pairing field. 

It is found that the LDA leads to a spatial variation of the gap near the surface 
of a nucleus, which is stronger than the one obtained in the HFB calculation. This 
is caused by the neglect of the proximity effects and the delocalized character of 
the single-particle wave functions close to the fermi energy. 

From the energy of the system, 

<E>=^n q E q , (5.35) 
i 

where n q = (1 + e Eq / T )~ 1 is the occupation number (at the temperature T) for 
the quasi-particle state q, one can obtain the specific heat of the system: 

C -i^. (5.36, 

This quantity turns out to be very sensitive to the pairing gap, and hence to the 
presence of the inhomogeneity induced by the presence of the nuclei in the crust. 
The overestimate of the gap, as obtained in LDA, leads to a specific heat of the 
system which is too large at low temperatures, as compared with the quantal 
result. Incidentally, a reliable estimate of the pairing gap is also important for 
the determination of the cooling time of the neutron star. 

In concluding this argument, it is worth mentioning that pairing correlations 
have been studied within a relativistic mean field approach based on a field 
theory of nucleons coupled to neutral (a and u>) and charged (p) mesons. The 
HF and pairing fields are calculated in a self-consistent way [36]. The energy gap 
is the result of strong cancellations between the scalar and vector components of 
the pairing field. It is found that the pair amplitude vanishes beyond a certain 
value of momentum of the paired nucleons; the estimated gap is in agreement 
with non-relativistic calculations of this quantity. 
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5.4 Neutron stars: massive quark matter 

While in the previous subsection the attention was focussed on the surface region 
of a neutron star, where the low density advocates nucleonic degrees of freedom, 
its high central density, up to 5 -f- 10 po, seems to imply the existence of a core 
of massive quark matter. The latter can be described in a variety of different 
models, among which we shall consider here the so-called Color Dielectric Model 
(CDM), which has been successfully employed in connection with the single 
nuclcon properties (structure functions, e.m. form factors) as well as for the 
nucleon-nucleon interaction. 

The CDM entails confinement of the quarks in ordinary hadronic matter 
but also allows to describe, in mean field approximation, a phase transition 
to deconfined quark matter[37, 38, 39]. The model is defined by the following 
Lagrangian: 



where x is the color dielectric field and U (a, tt) the "mexican hat" potential of 
the chiral sigma-modcl. This Lagrangian describes a system of interacting u, 
d and s quarks, pions, sigmas and a scalar-isoscalar chiral singlet field \. The 
latter is related to the fluctuations of the gluon condensate around its vacuum 
expectation value. The coupling constants are given by g u ^ = g(fn ± £3) and 
g s = (2fx — fn), where f n = 93 MeV and fx = 113 MeV are the pion and kaon 
decay constants, respectively, while £ 3 = f K ± — f K o = —0.75 MeV. Hence the 
model contains only two free parameters, g and A4, which arc fixed to the values 
g = 0.023 GcV, A4 = 1.7 GeV. Confinement is obtained via the effective quark 
masses, which diverge outside the nucleon. 

The CDM is employed to describe the deconfined quark matter phase in the 
neutron star, while the relativistic field theoretical model of Walecka is used for 
the hadronic phase. Applying Gibbs criteria to this composite system, Drago 
et a/.[38] find that the pure hadronic phase ends at 0.11 fm~ 3 while the mixed 
(quark and hadronic) phase extends up to 0.31 fm~ 3 . 

According to this calculation, neutron stars with masses in the range 1.3 < 
M/Mq < 1.54 and radii of about 9 Km arc found. A neutron star with total 
mass of 1.4 M Q will consist of a crust of pure hadronic matter, a <~ 1 Km thick 
region of mixed phase and a core of pure quark matter. 

In a recent work[39] the EOS of quark matter based on the CDM has been 
discussed, taking into account the effects associated with finite temperatures. 
Both the evolution of neutron stars and the onset of supernovae explosion cru- 
cially depend upon the EOS of matter at very high densities and temperatures. 
At finite T the EOS considered here shows a decrease of the pressure and of 




(5.37) 
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the adiabatic index, leading to a deconfmement transition for densities slightly 
larger than the one corresponding to nuclear matter saturation. The presence 
of a mixed phase region softens the EOS and could lead to a direct supernova 
explosion. At larger densities the EOS is stiff enough to support a neutron star 
compatible with observations. 



6 Quark matter 

Quantum Chromodynamics (QCD), the non-abelian theory of coloured quarks 
and gluons, is currently accepted as the theory of strong interactions, and its 
predictions have been tested in a variety of elementary particle reactions at large 
momentum transfers. Its behaviour in the high energy (or short distance) limit 
approaches the one of a non-interacting free field theory (asymptotic freedom), 
while at low energies (or large length scales) QCD is a non-perturbative field 
theory, quarks and gluons being permanently confined inside hadrons. 

At very high nuclear densities and/or temperatures it is believed that hadronic 
matter undergoes a phase transition to a deconfined phase of quarks and gluons, 
the so-called Quark Gluon Plasma (QGP). In this new phase hadrons dissolve, 
strong interactions become very weak and an ideal colour-conducting plasma 
of quarks and gluons is formed. In the QGP the long-range colour force is 
Debye screened due to collective effects and the quarks can only interact via 
a short-range effective potential. 

During the last decade the attention of nuclear and particle physics has been 
attracted by the possibility of producing this new state of matter in the lab- 
oratory, by means of ultra-relativistic heavy ion collisions. A large number of 
experiments have been carried out at the AGS (Brookhaven) and at the SPS 
(CERN). Both laboratories have planned considerable "upgrading" of the ex- 
isting facilities, the so-called RHIC (Relativistic Heavy Ion Collider), which is 
going into operation in 1999 and the LHC (Large Hadron Collider), whose pro- 
gram will be partly dedicated to high energy particle physics and partly to 
heavy ion collisions. The most recent measurements refer to S+U reactions at 
200 GeV/A and Pb+Pb at 158 GeV/A, reaching an (estimated) energy density 
of 2^-5 GcV/fm 3 . 

Many aspects of the transition from hadrons to deconfined quark matter, 
among which the order of the phase transition and the nature of experimentally 
observable signatures, are still under debate. Recently the NA50 experiment at 
CERN [40] has shown an anomalous suppression of the J/'P survival probability 
measured in Pb on Pb collisions [41], which might be ascribed to the formation of 
QGP at some stage of the collision. Other interesting signals could be revealed 
by the dilepton spectra, which are sensitive to the masses of vector mesons: the 
latter can be significantly altered by the presence of the QGP, an argument 
which is still widely discussed. 
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6.1 Quark deconfinement 

At high densities and temperatures, the interaction between quarks and gluons 
dresses their propagation, so that gluons develop an effective mass, which in 
turn produces screening of the long-range colour-electric forces. For example, 
the free gluon propagator D (u>,k) ~ (u 2 — fc 2 ) -1 is modified by summing an 
infinite chain of one loop insertions. One obtains [42] the following longitudinal 
gluon propagator 

D ^ = ¥I^F) (6 ' 38) 

Dl being a scalar function of the variables u) = fc° and k = |k| and el the so- 
called colour dielectric function. From explicit calculations, one can show that 
the static (w = 0) longitudinal colour fields are screened, being 



111 

k 2 e L (0,k) ~ fc 2 + g 2 T 2 = k 2 + X D 2 



°l(0, k) = -^-j—, = 17r —^ = t^^2 (6-39) 



which defines the Debye length Ad = (gT)^ 1 . 

The above mentioned features of the gluon propagator entail several impor- 
tant consequences, among which we focus here on the fact that the potential 
between two static colour charges (e.g. two heavy quarks) is screened in the 
quark-gluon plasma phase. Indeed the Fourier transform of Dl(0, k), eq.(6.39), 
yields the potential 

V QQ (r) ~ \e-rl^ (6.40) 

with screening length An — 0.4 fm at T = 250 MeV. This screening of the long 
range colour forces is believed to be at the origin of quark deconfinement in the 
high temperature phase. It also leads to the disappearence of the bound states 
of a charmed quark pairs (cc) in the QGP[43]. 

Fig. 10 illustrates the deconfinement effect of a thermodynamical enviroment 
of interacting light quarks and gluons on the interquark potential employed in 
ref.[44] to investigate the deconfinement conditions for heavy quarkonia.: 

V(r, T) = —pfn (l — ^ {T)r ) -^ iT)r (6-41) 

where fi{T) = l/An(T). The potential (6.41) reduces to the Cornell potential in 
the zero temperature limit (we remind that Ad — > oo for T — > 0); in addition 
to the one gluon exchange contribution (1/r) this potential also contains a phe- 
nomenological linear binding term, as suggested by lattice QCD calculations for 
the interaction between static (heavy) colour charges. 

It would be important to find some "QCD inspired" model, suitable to de- 
scribe both the confined (hadronic) and the deconfined (plasma) phases of a 
system of many quarks. This is obviously not trivial to achieve, due to the 
non-perturbative nature of the quark-quark interaction in the confined phase 
and also to the unavoidable approximations required to describe a many-body 
system, even when the interactions are known. 
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Fig. 10. The quark-quark potential, eq.(6.41) at T = 250° (dashed line) and T = 
(continuous line). 

In a recent work Alberico et al. [45] have developed a three-dimensional model 
for quark matter with a density dependent quark-quark (confining) potential, 
which allows to describe a sort of deconfinement transition as the system evolves 
from a low density assembly of bound structures to a high density free Fermi 
gas of quarks. Different confining potentials are considered, some of which suc- 
cessfully utilized in hadron spectroscopy, like the Cornell one. 

We find that a proper treatment of the many-body correlations induced by 
the medium is essential to disentangle the different nature of the two (hadronic 
and deconfined) phases of the system. For this purpose the ground state energy 
per particle and the pair correlation function are investigated. The latter can be 
obtained from the expectation value of the two-body density operator: 

9(r)= N{N p2 ~ 1} <^p 2 (|r 1 -r 2 l)i^>, (6.42) 

where {'P > is the exact (normalized) ground state of the system. For a system of 
strongly correlated pairs g(r) will show up a well localized peak at small values of 
r, while the Fermi gas pair correlation function has a fairly constant behaviour, 
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with the exception of small r values, where the Pauli principle prevents particles 
to be close to each other. 

As an example, Fig. 11 shows the pair correlation function derived from the 
Bethe-Goldstone wave functions of the screened Cornell potential: 

Vcomeii(r, p) = (-£ + br + K) e- c " r (6.43) 

where a, b, c and K are constants, while p is the density of the system. The 
free Fermi gas two-body correlation function is recovered for p ~ 0.6 fm -3 , in 
agreement with the value of the "transition density" which can be obtained from 
the corresponding equation of state. 




Fig. 11. The pair correlation function g(r) obtained with the potential (6.43) is 
displayed as a function of the relative interquark distance at various densities: 
p — 0.05 fm~ 3 (thick solid line), p — 0.1 fm~ 3 (dashed line), p = 0.2 fm~ 3 (dot-dashed 
line), p = 0.35 fm~ 3 (dashed line) and p — 0.6 fm~ 3 (dotted line). For the last density 
value the Fermi gas correlation function is also shown (thin solid line) . 



6.2 Finite size effects in the QGP 

Before concluding this Section it is worth mentioning that one of the major 
difficulties in the analysis of the relativistic heavy ion collisions stems from the 



Nuclear matter, nuclear and subnuclear degrees of freedom 



29 



numerous uncertainties in the thermostatistical evolution of the system. Among 
the various steps, the hadronization process is crucial in determining the final 
states measured in the detectors. 

Brink et al. [46] have studied the hadronization of a plasma, in order to estab- 
lish whether the final particles are evaporated from droplets of hadronic matter 
inside the plasma or rather from the whole volume of the plasma itself. The two 
situations can be distinguished on the basis of the entropy density of the system. 

Let us assume that at equilibrium (close to the critical temperature T ~ T c ) 
the volume of the plasma (Vtot) is filled by N droplets of QGP (each occupying a 
volume Vgiob) and a gas of relativistic pions. For a fixed total volume and energy 
density, the calculated entropy density of the systems turns out to be larger 
when the plasma fills a big, unique glob rather than several globs of smaller 
radius. This is partially illustrated in Fig. 12. Finite size corrections are taken 
into account in the density of levels which enters into the calculation of the 
entropy Fluctuations with respect to the number of globs are small, hence there 
is no tendency for a large glob to break spontaneously into smaller ones. 
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Fig. 12. Entropy density for a gas of quarks and gluons as a function of the glob radius 
R. The temperature is fixed to be T = 200 MeV. 



This argument concludes the tour which was foreseen in the introduction, 
showing how the various many-body techniques can be applied to strongly in- 
teracting systems, either of nucleons or of quarks, in an attempt to provide a 
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unified microscopic description of the intriguing, complex structure of nuclear 
systems. 

Acknowledgments 

I would like to thank Dr. A. de Pace and Prof. A. Molinari for enlightening 
discussions and to acknowledge their valuable help. 

References 

1. J.E. Amaro, M.B. Barbaro, JA. Caballero, T.W. Donnelly and A. Molinari, Nucl. 
Phys. A643, 349 (1998) 

2. R. Cenni, T.W. Donnelly and A. Molinari, Phys. Rev. C56, 276 (1997) 

3. W.M. Alberico, G. Chanfray, J. Delorme, M. Ericson and A. Molinari, Nucl. Phys. 
A634, 233 (1998) 

4. See the contribution of C. Ciofi degli Atti to these Proceedings 

5. More details can be found in the contribution of M.B. Barbaro to these Proceedings 

6. W.M. Alberico, M.B. Barbaro, S.M. Bilenky, J.A. Caballero, C. Giunti, C. Maieron, 
E. Moya de Guerra and J.M. Udias, Nucl. Phys. A623, 471 (1997) 

7. W.M. Alberico, M.B. Barbaro, S.M. Bilenky, J.A. Caballero, C. Giunti, C. Maieron, 
E. Moya de Guerra and J.M. Udias, Phys. Lett. B 438, 9 (1998) 

8. W.M. Alberico, M.B. Barbaro, S.M. Bilenky, J.A. Caballero, C. Giunti, C. Maieron, 
E. Moya de Guerra and J.M. Udias, e-Print hep-ph/9812388 (Dec. 1998) 

9. A. Fabrocini, Phys. Rev. C55, 338 (1997) 

10. R.B. Wiringa, V. Fiks and A. Fabrocini, Phys. Rev. C38, 1010 (1988) 

11. R. Schiavilla, V.R. Pandharipande and D.O. Riska, Phys. Rev. C40, 2294 (1989) 

12. Z. Meziani et al, Phys. Rev. Lett. 52, 2130 (1984) 

13. T.C. Yates et al, Phys. Lett. B312, 382 (1993) 

14. J. Jourdan, Nucl. Phys. A603, 117 (1996) 

15. J.E. Amaro, A.M. Lallena, G. Co and A. Fabrocini, Phys. Rev. C57, 3473 (1998) 

16. R. Cenni, F. Conte and P. Saracco, Nucl. Phys. A623, 391 (1997) 

17. T.S. Wahlhout, R. Cenni, A. Fabrocini and S. Fantoni, Phys. Rev. C54, 1622 
(1996) 

18. A. De Pace, Nucl. Phys. A (1998) 

19. M.B. Barbaro, A. De Pace, T.W. Donnelly and A. Molinari, Nucl. Phys. A598, 
503 (1996) 

20. M.B. Barbaro, R. Cenni, A. De Pace, T.W. Donnelly and A. Molinari, Nucl. Phys. 
A643, 137 (1998) 

21. B.D. Serot and J.D. Walecka, Int. J. Mod. Phys. E6, 515 (1997) 

22. P. Amorc, R. Cenni, T.W. Donnelly and A. Molinari, Nucl. Phys. A615, 353 (1997) 

23. A. De Pace, C. Garcia-Recio and E. Oset, Phys. Rev. C55, 1394 (1997) 

24. J.A. Oiler, E. Oset, A. De Pace and P. Fernandez de Cordoba, Phys. Rev. C57, 
1404 (1998) 

25. A. De Pace, in Proceed, of the Int. Conf. on Hypernuclear and Strange Particle 
Physics, Brookhaven Nat. Lab., Oct. 1997 

26. CM. Kormanyos et al, Phys. Rev. C51, 669 (1995) 

27. G. Garbarino, W.M. Alberico, A. De Pace and A. Ramos, in preparation 



Nuclear matter, nuclear and subnuclear degrees of freedom 



31 



28. G. Co and A. Fabrocini, work in progress 

29. V. Baran, M. Colonna, M. Di Toro, A.B. Larionov, Nucl. Phys. A632, 287 (1998) 

30. M. Di Toro, M. Colonna, V. Baran, A.B. Larionov, Int. Conf. on Giant Resonances, 
Varenna 1998, Nucl. Phys. A; see also the contribution of M. Di Toro et al. in these 
Proceedings 

31. Madappa Proakash, I. Bombaci, Manju Prakash, P.J. Ellis, J.M. Lattimer and R, 
Knorren, Phys. Reports 280, 1 (1997) 

32. M. Baldo, I. Bombaci and G.F. Burgio, Astron. Astrophys. 328, 274 (1997) 

33. I. Bombaci, Phys. Rev. C55, 1587 (1997) 

34. F. Barranco, H. Esbensen, E. Vigezzi and R.A. Broglia, Phys. Lett. B390, 13 
(1997) 

35. F. Barranco, R.A. Broglia, H. Esbensen and E. Vigezzi, e-print nucl-th/9806056 

36. F. Matera, G. Fabbri and A. Dellafiore, Phys. Rev. C57 (1997) 

37. A. Drago, M. Hjorth- Jensen and U. Tambini, Prog. Part. Nucl. Phys. 36, 407 
(1996) 

38. A. Drago, U. Tambini and M. Hjorth-Jensen, Phys. Lett. B380, 13 (1996) 

39. A. Drago and U. Tambini, preprint, March 1998 

40. M. Gonin, NA50 Collaboration, Nucl. Phys. A610, 404c 

41. D. Kharzeev, C. Lourenco, M. Nardi and H. Satz, Z. Physik C 74, 307 (1997) 

42. J. Kapusta, Nucl. Phys. B148, 461 (1979) 

43. T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986) 

44. F. Karsch, M.T. Mehr and H. Satz, Z. Physik C 37, 617 (1988) 

45. W.M. Alberico, P. Czerski and M. Nardi, Eur. Journ. Phys. A (1999), in press 

46. D.M. Brink and Lo Monaco, J. Phys. G 24, 867 (1988) 



